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Abstract 

I explicitly work out closed form solutions for the optimal hedging 
strategies (in the sense of Bouchaud and Sornette) in the case of Eu- 
ropean call options, where the underlying is modeled by (unbiased) hd 
additive returns with Student-t distributions. The results may serve 
as illustrative examples for option pricing in the presence of fat tails. 



1 Introduction 



Black and Scholes option pricing theory relies on the existence of a dynamic 
trading strategy in the underlying asset that exactly replicates the option 
holder's claim [1-4]. There is no risk for the writer, and the fair option price 
can be interpreted as the cost of the perfect hedging strategy. 

Real world markets are, however, not complete, i.e., there are no riskless 
hedging strategies that would provide unique option prices. The incomplete- 
ness stems from a number of reasons. E.g., stock market returns are not 
log- normally distributed (as assumed in the Black Scholes analysis). It is 
widely agreed that marginal return distributions have fat tails, i.e., large 
events occur with significantly higher frequency than to be expected from 
the log-normal assumption. 

An interesting approach to option theory for incomplete markets has 
been worked out by Miiller ||, Follmer and Sondermann ||, Schal and 
Schweizer ||, and also by Bouchaud and coworkers ||, [UJ, see also [11-15]. 
At the heart of these approaches lies a minimization of the risk inherent in 
writing the option. Different risk measures have been suggested, including 
the global and local variance of the cost process [5-7] and the variance of the 
global operator wealth 11. 

In this article, I will explicitly work out closed form solutions for the 
optimal hedging strategies (in the sense of Bouchaud and Sornette) in the 
case of European call options, where the underlying is modeled by (unbiased) 
iid additive returns with Student-t distributions. The results may serve as 
illustrative examples for option pricing in the presence of fat tails. Further- 
more, since Student distributions have been considered as acceptable models 
of market returns [16-20], the results presented here could also be useful for 
practical purposes. 

This article is organized as follows: In section 2, I setup the model and 
notation. A number of useful formulae is derived that are of significant help 
in the explicit calculations of later sections. In section 3 we reproduce a few 
(well known) results for Gaussian returns. In section 4 I derive the closed 
form solutions^ for the Student-t case. Illustrative example results are given 
in section 5. Conclusions follow. Some technical parts are deferred to three 
appendices. 

1 Closed form here means that the result is presented as a finite sum of terms built from 
standard functions. 



1 



2 Model Setup and Basic Formulae 



Consider a European call option, with maturity at T = Nt, and strike 
price x s . Trading and hedging is only possible at discrete times t — kr, 
k = 0,...,N. The underlying price process is given by x k = x(kr). At 
expiry the option holder's claim is given by 

/(a^jv) = max(xjv — x s , 0) . 

For the sake of simplicity we assume that the (risk free) interest rate is zero. 

After writing an option, the bank will perform the following (hedging) 
actions, k — 0, . . . , N—l: Trade the underlying such that the bank's portfolio 
contains 4> k pieces of stock, k = N: Sell the portfolio and satisfy the option 
holder's claim. 

Following Bouchaud and Sornette |§, we write down the global cost bal- 
ance for the bank (per option issued) :Q 

N-l 

C = max(xAr - x s , 0) - 4>k(x k ) (x k+1 - x k ) . (1) 

fc=0 

If it were predictable (deterministic), C would be the correct price of the 
option (ignoring transaction action costs etc.). However, as a function of the 
underlying, C is a stochastic variable. The fluctuations of C will in general 
depend on the choice of the hedging strategy <p. In the Bouchaud Sornette 
approach, on chooses in order to make C as determininistic as possible, 
e.g., by demanding that the variance of C is minimal. The optimal hedging 
strategy 0* is then determined by 

# 2 (0*) = (C 2 ) - (C) 2 = minimal , (2) 

where (.) denotes averaging over the price process. For sufficiently small 
residual risk, e.g., for Rj (C) <K 1, it is reasonable for the bank to define the 
option price as 

price = (C) + premium for the residual risk . 

The extra premium for the residual risk will of course depend on the bank's 
risk preferences. 

2 We have not included the option premium in the balance. 
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Let us now assume the following (simplified) model for the returns of the 
underlying asset: 

x k+ i - x k = r k+1 , k = 0, . . . , N - 1 , 

where the r's are iid with probability density function (pdf) p(r). We require 
that p should have a finite variance, defined by0 



a 2 = J drp(r) r A 



However, we are not assuming that p is Gaussian or log-normal. The stock 
price x k can be written as 

k 

x k = x + ^2r k = x + R k . 
i=i 

The probability measure is defined by 

(F) = / dr x p(n) . . . / dr N p(r N ) F(r h ...,r N ). 



We will assume that p(r) is even, i.e., p(r) = p(—r). Then the expectation 
value of C is unaffected by the hedging strategy 0: (C) = (/). In the present 
context it is given by 

/OO /"OO 
dt P N (t) max(x + t - x s ,0) = / dtP N (t)(( + t) 
-oo J — £ 

We have defined 

C = Xq — XS . 

At k = the option is in-the-money, at-the-money, or out-of-the-money, for 
( > 0, ( = 0, or ( < 0, respectively. 

is the pdf of Ri, i.e., the /-fold convolution of p. It can be conveniently 
accessed via Fourier transformation methods. Introducing the characteristic 
function p(q) through 

p(q) = dr e iqr p{r) , 



we have 

PM) = P(q) 1 ■ 

integrals without explicit limits should always be extended over the whole real axis. 
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It is not difficult to show that the hedging strategy minimizing the vari- 
ance of C can be expressed as 

4>l{ x k) = -J / drp(r)r J dt Pjv-fc-i(t) f(x k + r + t) . 

It is demonstrated in Appendix A that 

1 f dq sin(g£) 1 d , N _ k 



where <fi* k = \ + V'fc? an d £ = x k — x a . 
The residual risk can be expressed as 



JV-l 

2/ j* 



R z {<t>*) = (f) - (f) 2 - ° 2 E / - *o) 4>l(xf • (4) 

fc=0 



Here, (f 2 ) — (f) 2 is the unhedged risk. (/) and (f 2 ) can be computed with 
the help of the following relations, derived in Appendix B: 



(/> = C (§ + MC)) + ^(C) , 

(f 2 ) =C 2 (f + Li(C))+C£ 2 (C)+(^a 2 + L 3 (C) 



(5) 



where 



Li(c) =J ^sinM~ (g) iV 

L 2 (C) =-/gcos(Cg)^%r, (6) 



q dq 

The second term in eq. (|3j) can be rewritten as 

jV-1 



iV-l ,1 

fc=0 7 v / 



In the examples to be discussed below we shall evaluate this term by per- 
forming the ^-integration numerically, starting from closed form expressions 
for Pk and t/>£. 
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3 The Gaussian Case 



As a first application of the formulae derived in the previous section, and 
for later comparison with the Student-t case, let us consider the case of 
Gaussian returns, i.e., p(r) = 1 j^l-no 2 e~ r ^ 2<T \ The characteristic function 
is p(q) = e~2 a q . The integrals for and Lj are elementary. With of = I a 2 
we find 

r k>G = ^rf (e/ y/toj^) . (7) 

For 4>l G = \ + i ) kG one ^us obtains 4>* kG = N(£/aN-k), where N denotes 
the cumulative distribution function of the normal distribution. For the Li 

we get 



L hG =± erf CAM 



C 2 /(2<t2 ) 



J^2,G — &N 



L 3 , G =la%erf((/J2a 



For a complete discussion of the Gaussian case see [0, 10 



4 Options with Student-t Underlying 

We shall now consider the case that the return distribution is of the Student-t 

type, 




a > is a scale parameter that we shall later use to fix the standard deviation. 
The latter exists for p > 2. For p an odd integer, one can derive relatively 
simple expressions for the characteristic functions p^> , thus making feasible 
closed form solutions for the hedge functions (f>%. Let us therefore concentrate 
on the case p = 2m + 1, with m = 1, 2, . . .. By an elementary integration 
one can show that 



A straightforward but somewhat lengthy calculation (employing the theorem 
of residues) yields 

p( 2m+1 \q) =e~ aM T [m] (a\q\), 



where 

m / ™ \ / / o™ \ to „\k 

T [m \q) = Y. 



in \ I I 2m \ (2qY 



k =o v k J I \ k ) k\ 
is a polynomial of degree m. Here are a few examples: 

TW(q) =1, 
TW(g) =l + q, 
T®(q) =l + q+\q\ 

In Appendix C eq. (|3]) is used to derive a closed form expression for ijj 
The result is 

1 H m(N-k)-l A [m,N-k] _ -.x , sin A arctan f^)") 

= - arctan -A- + £ ^ — ^ 

Here we have defined 



7/i * 



a a 



The coefficients Aj' N fc ' are defined by 

m(N-k)-l 

T [m] ^N-k-1 T [m-1] (g) _ £ ^KiV-fc] ?j 

i=o 

Calculations very similar to those presented in Appendix C help us to find 
explicit expressions for the Lj, defined in eq. (^). 
For Li we get 

L [m] = r dq sin(C'g) c -N„ T \m], sN 
1 Jo 7r g 

with £' = C/a. Expanding 

mff 

T lm]( q) N = J2 B [m,N] ^ 
1=0 
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we find 

L M = I arctan C f (/-!)! sin (/ arctan (^)) 

1 7T N 71 (('2 + jV 2 ) Z/2 

Similarly, for L2 one obtains 

L2 = E^L r^L CO s(Cq)e- Nq T^(q) N ~ 1 T^- 1 \q). 
a Jo ix 

Using that 

f°° b ^(0) ( T 1 

/ rfx e - ^ cos(7x) = — — — — cos a arctan — 

Jo {(3 2 + 7 2 ) Q/2 V /V 

one arrives at 

w _ Na 2 A [ ™> N] l\ cos ((/ + 1) arctan (£)) 

The calculation of L 3 goes along the same lines. We find 

/l C / ^i ^MVl (/_!), sin (/ arctan 

Ln = Na — arctan — + > j^r- 

3 [ir N n (C' 2 + N 2 f 2 

For the calculation of the residual risk we shall also need an expression 
Pk, the fc-fold convolution of p. The result is 

mfc B \mM ^ cog fn + j\ arctan (*)) 
a P k {at) = V ^ ^ ' n^,} JJ ■ 



for 



5 Example Results 

In this section we will present some example results for the Student-t distri- 
bution, in comparison with the Gaussian case. Let us start with a look at 
the hedge functions for different values of m. We choose N = 7 and k = 0. 
If you think of days as the hedging and trading intervals, this corresponds 
to one week before expiry. Figure |l] shows the results for m = 1, 2, and 3, 
together with the Gaussian results. i/jq is plotted as function of (/a, where 
a denotes the standard deviation (per day). 
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Figure 1: The N = 7, k = hedging function ipQ for Student-t 
underlying with m = 1, 2, and 3 (bottom to top). The uppermost 
curve is the hedge function for the Gaussian case. 



With increasing m, the curves converge quickly towards the -^-function 
of the Gaussian returns. Closer to expiry the difference to the Gaussian 
behaviour is more significant, as demonstrated in Figure 0, which shows ipo 
for the same set of m-values, however with N — k = 1, i.e., one day before 
maturity. 

Let us next consider the behaviour of (C). Figure ^| shows this quantity 
for N = 7, for m — 3, and for the Gaussian, as a function of C/cr. Around 
C = 0, the Student result is a little bit below the Gaussian, farther out it is 
a very little bit above it ("volatility smile"). 

Figure |] shows again results for N = 7, (C) as a function of C/cr, together 
with the risk (C 2 ) — (C) 2 , plotted with the help of error bars. The upper plot 
gives these quantities for the unhedged case, i.e., for = 0. For the lower plot 
the optimal hedging was performed. In both cases, the data with the smaller 
risk belong to the Gaussian, the other ones to the Student-t with m = 1. 
Note that there remains quite some risk even after the optimal hedging. It 
is only for the Gaussian that this residual risk can be removed in the limit 
of continous hedging. 
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Figure 2: The TV = 1, k = hedging function V'o f° r Student-t 
underlying with m = 1, 2, and 3 (bottom to top). The uppermost 
curve is the hedge function for the Gaussian case. 
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Figure 3: (C) for N = 7, Student-t, with m = 3, in comparison with 
the Gaussian. The Student-t result is below the Gaussian for small (. 
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Figure 4: For N = 7: (C), together with (C 2 ) - (C) 2 , plotted with 
the help of error bars. Top: unhedged, bottom: optimally hedged. 
Bars: Gaussian, crosses: Student-t with m = 1. 
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Figure 5: Upper curve: residual risk for N = 7, Student-t with 
m = 1, hedged with Gaussian hedge function. The lower curve is the 
residual risk in the case of optimal hedging. 



It is interesting to visualize the "hedging error" occuring when one hedges 
the Student-t underlying with the Gaussian hedge functions given by eq. (^)). 
The upper curve in Figure [| denotes the residual risk for N = 7, Student-t 
returns with m = 1, hedged with the Gaussian hedge function (corresponding 
to the same standard deviation). The lower curve is the residual risk in the 
case of optimal hedging. The difference is not dramatic, but appreciable. 

For non-Gaussian distributions, the variance is, of course, not always an 
appropriate measure of risk. It is therefore instructive to look at the full 
pdf of the quantity C. Given the hedge functions, it can be easily computed 
by Monte Carlo. In Figure || we see histograms (from 1 million samples of 
the whole hedging process) of C (N = 7) for the Gaussian and the m — 1 
Student-t returns, both optimally hedged. In the lower part of the figure 
(with log scale), we show additionally the function l/(C/a) 4 , demonstrating 
that the power tail of the Student-t distribution is still around. For compari- 
son: The amplitude in front of the 1/x 4 term in the 7-fold convolution of the 
m = 1 Student-t (with standard deviation one) is 0.24062. If one assumes 
that the option writer would sell the option for (C) = 0.96, his value at risk 
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Figure 6: The pdf's of C, for N = 7. Top: linear scales, bottom: 
vertical log scale. The distributions with the fatter tails belong to 
Student-t with m = 1, the remaining stems from Gaussian returns. 
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(VaR(5 %)) would be 0.98, i.e., more than 100 percent of the option price! 

6 Conclusions 

In this paper, closed form solutions were derived for the optimal hedge func- 
tions in the case of unbiased Student-t returns. I hope that the results prove 
useful as illustrative examples for option pricing in the presence of fat tails, 
and as starting points for further investigations. 
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Appendix A: Representations of <jf k 

We start from 

<f>t(x) = ^J drp{r)r J dt P N - k -!(t) f(x + r + t) , 
with f(x) = max(i-a; s , 0). The claim function can be represented as follows: 

dq e iq ( x ~ Xs ^ 



/(*) = - / 



2% (q — ie) 2 

To see this one evaluates the integral using the theory of residues. 



r p iqX d 



= -2tt 0(A) A = -2ix max(A, 0) . 



q=ie 



Here, 9(\) denotes the Heaviside step function. With xu — x s = £ we thus 
have 

dq e iq ^ +r+t) 
{q-ief 

Exchanging the order of integration, we obtain 



1 r r r An e n\K+r+t) 



The r-integral can be written as \j-p{q)-, whereas the t- integral is the same 
as p(q) N ' k ^ 1 . Putting things together we obtain 



6* = 1 [ ^ ^ ± f,(a) N 1 

9k (N -k)a 2 J 2m (q - ie) 2 dq m> 



(N — k) a 2 J 2m (q-ie) 2 dq 
To proceed further, we write e iq ^ = cos(g£) + i sin(g£). Then 
/ * 1 f d( l cosg^ d N „ k 

^ = -W^W^J^¥^Wdq p(q) () 
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with 



% ~ (N-k)a^J 2n q q dq P[q) 



In the expression for ip* we have already performed the limit e — > 0. The first 
term of eq. (H) yields 1/2. To see this we observe that the term is independent 
of £: Differentiating with respect to £ makes the integrand an odd function 
of q and the integral vanish. Without changing the integral we can therefore 
let £ = 0. Furthermore 



(N-k) a 2 J 2ni (q-ie) 2 dq rKH) 



(N-k) a 2 I dxP N - k (x) X J 2 l ( q } ie) 2 
(N-k) a 2 io ^ ^W-fc(^) ^ 2 



This proves the relation 0£ = ~ + r/>£. 

Appendix B: Computation of (f n ) 

The equation 

(r> = Jdtp N (t)6(c+t) (c+tr 

can be rewritten with the help of the following integral representation of the 
Heaviside step function: 

0{x)= fdQ ^ 



2-ni q — it 

where e > is to be considered infinitesimal. We obtain 

dq 1 



2ni q — it 

By elementary operations this can be rewritten as 

dq 1 fid 
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Working out the derivatives for n = 1,2, we obtain the representation given 
in eqs. (|5|,|6|), provided we recognize the identities 



and perform an integration by parts in one of the integrals. Eqs. ([10| , |TT|) 
can be proved by first observing that both expressions do not depend on 
(. Letting ( = thus leaves the integrals invariant. Transforming back to 
t-space we recognize integrals of P/v(*) an d -P/v(*) & \ respectively, for t from 
to oo. 



Appendix C: Computing tp* for Student-t 



It is easy to derive (by partial integration) the following recursion relation: 




(10) 




(11) 



We start from eq. @ in the form 





Employing this recurrence and eq. (|5|), we obtain 




Performing the substitution q — > q/a and defining 



a 



€ x k -x. 



s 



a 



one obtains 




(12) 
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The product of the two T-factors is a polynomial of degree m (N — k) — 1: 

m{N-k)-l 

T [m] T^-V (q) = A [ ™' N ~ k] q l . 

1=0 

E.g., for m — 1, the ^-coefficients are given by 

For larger m the expressions become more involved. They can however easily 
be obtained with the help of a computer algebra program. 
Eq. (|I2|) becomes 



-i m(7V-fc)-l „ ■ ( H\ 

7T Jo g 



-(N-k)q q l 



1=0 

This integral can be solved in closed form (Gradstein/Ryshik, 3.944): 

7(a, 7) = y o ^x Q_1 e~ l3x sin(7x) = 2 + 2 sin la arctan -J . 

7 (13) 
We have to identify a = I, j3 = N — k, and 7 = Observing that 

7 

lim /(a, /3, 7) = arctan — , 

we finally arrive at 

1 £' m ^_ fc )- 1 aJ™^-" 1 (Z - 1)! sin (l arctan 
Wk =~ arctan l\r T +2^ 



7T N-k ^ it (£« + (jv - A;) 2 ) i/2 

Note that this can alternatively be represented as 



1 £' A! miJV " fel (Z-l)! Im(JV-ife + iO 



fc vr N-k £i 7T (£'2 + ( N _ k yy 

One might wish to expand the imaginary part 

[(/-l)/2] 



[l«-J-J/2j / , \ 

Im(a + z6)'= ]T ( ■ ) (-I)' a l - 2 ^H 2j 
j=o \ 3 J 



+1 



thus demonstrating that the corrections to the leading arctan behaviour de- 
cay like for large [7] denotes the largest integer smaller or equal to 
3- 
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